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Abstract 

Theoretical  and  numerical  studies  of  metamaterials  are  addressing  the  question  of  whether  true 
optical  magnetism  is  possible.  Lattice-based  structures  that  have  symmetries  leading  to  unique 
descriptions  of  homogenized  permeability  are  being  evaluated.  Current  work  on  conduction  current 
magnetism  in  carbon  nanotubes  is  being  extended  to  more  general  materials  with  consideration 
of  displacement  current  effects.  Numerical  results  are  being  viewed  both  in  terms  of  spatial  dis¬ 
persion,  associated  with  the  structure  of  the  material,  and  equivalent  circuit  concepts.  Promising 
results  from  this  work  will  guide  plans  for  an  experimental  demonstration  of  optical  magnetism. 
Establishing  that  optical  magnetism  is  possible  would  be  of  enormous  importance,  both  in  terms 
of  the  fundamental  science,  and  also  the  myriad  of  applications  that  could  be  revolutionized  by  a 
magnetic  metamaterial. 
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I.  INTRODUCTION 


In  this  report,  we  present  the  case  for  homogenized  magnetism  from  arrays  of  carbon 
nanotube  coils.  There  are  existing  fabrication  methods  for  such  coils,  making  the  issue  of 
the  frequency  range  for  significant  permeability  important  and  timely.  We  present  a  study 
of  homogenization  of  LiNbC>3  (LN)  and  metallic  structures  to  obtain  the  effective  consti¬ 
tutive  parameters  from  inversion  of  scattering  parameters  S2 1  and  Su  [13].  The  structures 
consist  of  periodically  arranged  lattice  cells  of  dielectric  or  metallic  patterns,  with  lattice 
constant  on  the  order  of  100  /irn  with  a  view  to  evaluation  at  terahertz  frequencies  in  a 
scaled  picture  for  optical  magnetism.  Scattering  parameters  (Su  and  S21)  were  obtained 
from  finite-element  (FEM)  simulations  of  the  structures  with  an  incident  monochromatic 
plane  wave.  A  procedure  for  inverting  these  scattering  parameters  based  on  the  formulation 
of  Nicolson  and  Weir  was  applied  to  extract  the  homogenized  constitutive  parameters.  Over 
the  frequency  involved  in  the  study,  the  effective  wavelength  in  the  structure  remained 
larger  than  about  10  times  the  crystal  periodicity,  which  justifies  the  homogenization.  The 
results  indicate  that  both  paramagnetic  and  diamagnetic  responses  can  be  achieved  through 
structured  patterning  of  LN  as  well  as  in  metallic  structures.  We  took  into  consideration 
the  fact  that  most  metamaterials  are  anisotropic  and  exhibit  chirality,  and  investigated 
the  dispersion  relations  of  plane  wave  propagation  in  chiral  (bi-isotropic)  and  the  more 
general  bianisotropic  media.  The  plan  for  subsequent  studies  includes  an  extension  of  the 
S-parameter  homogenization  method  to  allow  the  extraction  of  the  full  tensor  components 
of  er,  Jirl  and  in  addition  the  chirality  parameter  tensor  k. 


II.  MAGNETISM  FROM  CARBON  NANOTUBE  COIL  ARRAYS 

Our  preliminary  studies  suggest  that  THz  and  perhaps  optical  magnetism  is  possible  with 
carbon  nanocoil  mixtures  which  have  been  fabricated,  and  that  we  may  be  able  to  achieve  a 
resonance  in  the  infrared  spectral  regime.  Our  results  that  are  being  prepared  for  submission 
to  a  journal  are  summarized  here. 

Carbon  nanotubes  (CNTs)  and  nanowires  (CNWs)  have  been  the  subject  of  intense 
research  in  the  last  couple  of  decades  because  of  their  unusual  mechanical  and  electrical 
properties.  In  the  last  decade,  interest  in  more  exotic  nanostructures  such  as  CNW  rings  and 
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coils  has  evolved  both  due  to  academic  curiosity  and  to  the  great  potential  such  structures 
hold  in  electrical  and  mechanical  applications  [2,  4,  6,  14,  25].  However,  while  it  has  been 
known  for  some  time  that  growing  nanotubes  and  nanowires  in  increased  ambient  entropy 
(achieved  by  ultrasonic  and  microwave  chemical- vapor  deposition  (CVD),  or  injection  of 
sublimated  compounds)  results  in  coiling  [9,  10],  the  necessary  understanding  and  technology 
has  only  recently  matured  to  a  level  where  such  nanostructures  can  be  synthesized  [1].  This 
understanding  stems  from  a  unified  thermodynamic  model  that  explains  the  coiling  using  the 
energy  minimization  principle.  This  model,  together  with  recent  technological  improvements 
[11,  23,  24]  in  the  CVD  synthesis  of  coils,  paves  the  way  to  achieving  coils  with  high  yield  and 
high  precision  at  increasingly  small  coil  dimensions  (see  Fig.  1).  Those  achievements  lead  us 
to  believe  that  a  volume  containing  small  enough  nanocoils  at  a  high  enough  density  for  use 
as  an  effective  medium  may  produce  negative  permeability  at  THz  and  perhaps  some  optical 
frequency  ranges,  and  that  this  could  be  manufactured  using  today’s  CVD-based  methods. 


(a)  (b)  (c) 

FIG.  1.  (a)  [I  Sz  II]  Low  and  high  magnification  scanning  electron  microscopy  (SEM)  pictures 

showing  an  as-grown  forest  of  coil-shaped  carbon  nanotubes  (CCNT)  synthesized  by  use  of  thermal 
CVD  with  a  bi-metal  catalyst  formed  from  a  mixture  of  ferrocene  and  indium  isopropoxide  as  coiling 
seed.  [Ill  Sz  IV]  Low  and  high  magnification  transmission  electron  microscopy  (TEM)  images  of  a 
loose  single  bundle  of  coil-shaped  carbon  nanotubes  (from  [3]).  (b)  [I  Sz  II]  SEM  pictures  of  coiled 
features  of  CCNTs  synthesized  via  CVD  from  pyridine  and  sublimated  melamine  injection  (from 
[11]).  (c)  TEM  picture  of  carbon  CCNTs  synthesized  via  microwave  CVD  with  fine  Ni  particles  as 
catalyst  (from  [22]). 

Before  studying  the  effective  electromagnetic  parameters  of  a  volume  containing  a  multi- 
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tude  of  nanocoils,  it  is  important  to  pay  attention  to  the  subtle  issue  of  chirality.  Most  meta- 
materials  exhibit  only  simple  constitutive  relations  of  the  form  D  =  e^e-  E  and  B  =  /i0/i  •  H. 
where  E  and  H  are  the  electric  and  magnetic  fields,  respectively,  D  and  B  are  the  electric 
and  magnetic  flux  densities,  respectively,  eo  and  yU0  are  the  free-spaee  permittivity  and  per¬ 
meability,  respectively,  and  e  and  Ji  are  the  effective  relative  permittivity  and  permeability 
dyadics  of  the  metamaterial,  respectively.  However,  coils  are  asymmetric  with  respect  to 
mirror  reflection,  which  means  that  a  substance  comprising  coils  could  be  chiral,  and  hence 
exhibit  more  complicated  bi- anisotropic  constitutive  relations  of  the  form 

D  =  e0e  E  —  jy/JL 0e0  RT  •  H  B  =  jy/fj,0eo  re  ■  E  +  /i0/i  •  H,  (1) 

where  R  is  the  chirality  parameter  dyadic  relating  the  magnetic  held  to  the  electric  flux  and 
vise  versa  [18].  While  the  design  and  application  of  chiral  materials  has  been  thoroughly 
studied  during  the  last  decade  of  the  20th  century,  this  work  had  the  goal  of  attaining 
chirality  [17].  On  the  other  hand,  we  are  interested  in  suppressing  chirality  while  retaining 
the  large  local  magnetic  effect  of  each  coil.  By  doing  so,  it  is  possible  to  obtain  a  simple 
effective  medium  whose  quasi-static  resonance  will  have  a  large  enough  influence  on  the 
effective  permeability  to  drive  its  values  into  the  negative  range.  The  quasi-static  resonance 
is  due  to  the  inductance  of  the  coil  and  the  capacitance  between  coil  turns.  In  order  to 
suppress  chirality,  we  suggest  a  cellular  structure  where  each  cell  comprises  a  single  coil, 
and  each  pair  of  neighboring  cells  are  of  different  handedness  (i.e.,  left-hand  and  right-hand 
coils,  as  in  Fig.  2(a)).  This  structure  will  be  on  average  symmetric  with  respect  to  mirror 
reflection,  and  thus  non-chiral.  To  obtain  an  initial  estimate  of  the  effective  properties  of 
such  a  structure,  we  studied  the  homogenized  permeability  of  a  super  cell  of  8  CNW  coils, 
shown  in  Fig.  2(b),  in  a  doubly  periodic  lattice.  We  used  a  vector  electromagnetic  solution 
and  a  parameter  extraction  procedure  similar  to  the  one  used  in  [19],  as  a  function  of 
material  thickness,  in  order  to  establish  converged  bulk  parameters.  To  simplify  our  initial 
calculations  and  obtain  results  that  are  dependent  only  on  the  ratio  between  dimensions 
and  wave-length,  we  used  classical  graphite  properties  for  the  coils  and  assumed  a  small  skin 
depth  compared  with  coil  wire  thickness.  This  simplification  allowed  us  to  concentrate  on  the 
choice  of  coil  and  cell  dimensions  that  will  decrease  the  resonant  frequency  while  keeping  the 
coil  to  coil  interaction  weak,  so  that  the  effective  medium  parameters  will  converge  even  for 
small  material  thickness.  The  result  in  Fig.  2(c)  shows  the  real  and  imaginary  components 
of  the  homogenized  permeability  as  having  a  resonance,  and  that  it  is  a  metamaterial.  It 


4 


would  therefore  appear  that  a  mixture  of  CNT  coils  could  provide  negative  ji  at  THz  and 
perhaps  optical  frequencies  with  an  appropriate  geometry.  We  also  found  that  the  resonant 
frequency  in  the  homogenized  result  is  insensitive  to  source  polarization  and  orientation, 
which  is  very  important  for  random  mixtures. 


Right  handed  coil  Left  handed  coil 


a  a 


a  a 


(a)  (b)  (c) 

FIG.  2.  (a)  Left-handed  and  right-handed  carbon  coils  that  are  placed  in  a  x  a  x  a  cells.  The 

coil  diameter  and  length  is  0.5a,  and  the  wire  diameter  is  0.02a.  (b)  Non-chiral  unit-cell  of  peri¬ 
odic  metamaterial  that  comprises  4  left-handed  coils  and  4  right-handed  coils,  (c)  Real  (//)  and 
imaginary  (pw)  components  of  the  effective  permeability  of  the  coil  metamaterial,  obtained  using 
a  scattering  model. 


A  more  accurate  approximation  of  the  performance  of  CNT  coils  is  obtained  by  mod¬ 
eling  the  tubes  with  a  surface  conductivity  given  by  the  Landauer  formula  for  a  quantum 
conductor 


Qs,  tube 


GqMLq 
2t tR 


(2) 


where  Go  is  the  conductance  quantum,  M  is  the  number  of  modes  carried  by  the  wire, 
Lq  is  the  mean  free  path  length  in  the  wire,  and  R  is  the  tube  radius.  Based  on  the 
scaling  principle,  it  can  be  shown  that  the  resonance  frequency  is  independent  of  the  surface 
conductivity  value,  and  the  CNT  coil  effective  permeability  can  be  accurately  predicted  by 
a  simple  Lorentzian  resonance  model  in  the  angular  frequency  oj, 

A/zcn2 


/U  —  1  + 


cup  —  cu2  +  j2'ycu: 


(3) 


that  is  fitted  for  the  dipole  strength  A/i,  the  resonance  angular  frequency  u>o,  and  the  damp¬ 
ing  coefficient  7.  Figure  3  shows  the  effective  permeability  obtained  from  full-wave  simula¬ 
tions  and  by  the  model  in  (3).  Within  the  regime  of  weak  interaction  between  coil  elements, 
the  coefficients  of  the  model  in  (3)  have  a  simple  dependency  on  geometrical  properties  of 
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FIG.  3.  Real  (//)  and  imaginary  (fi")  components  of  the  effective  permeability  of  a  coil  meta¬ 
material  based  on  an  improved  tube  model  for  coils  with  different  number  of  turns.  The  dashed 
plots  are  obtained  using  scattering  parameters  extracted  from  full-wave  simulations  and  the  solid 
plots  are  obtained  from  the  resonance  model  in  (3) .  The  unit-cell  geometry  used  in  the  simulations 
consists  of  a  50  nm  diameter  coil  with  a  16  nm  pitch  that  is  set  in  the  middle  of  a  unit-cell  box 
that  is  constantly  kept  1.5  times  the  size  of  the  coil  in  all  dimensions.  The  surface  conductivity  of 
the  coil  is  based  on  the  Landauer  formula  in  (2)  with  M  =  2  modes,  a  mean  free  path  length  of 
Lq  =  150/xm,  and  tube  radius  of  R  =  0.5nm 

the  tube,  and  can  give  invaluable  insights  into  the  physics  governing  the  CNT  coil-based 
metamaterial  performance.  The  model  can  also  be  used  to  decouple  the  design  process  of 
CNT  coil-based  metamaterials  from  the  CNT  conductance  properties  that,  for  frequencies 
above  low  THz,  are  still  under  active  scientific  investigation  in  the  carbon  nanotube  research 
community. 

An  aspect  of  importance  for  the  applicability  of  nanocoil  based  metamaterials  is  their  loss 
tangent.  While  all  resonant  based  metamaterials  exhibit  some  degree  of  loss,  our  calculations 
show  that  this  loss  is  not  necessarily  significant  in  the  frequency  in  which  the  effective 
permeability  of  a  nanocoil  based  medium  is  negative.  It  is  interesting  to  note  that  while  some 
have  shown  that  substances  comprising  carbon  nanotubes  can  behave  as  a  dark  material, 
it  is  mainly  due  to  a  specific  choice  of  effective  medium  parameters  so  as  to  suppress  the 
substance’s  reflectance,  while  the  effective  absorption  length  of  such  substances  is  actually 
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very  long  [5]. 

Figure  4  shows  the  dependency  of  the  extracted  resonance  parameters  on  the  number  of 
modes  in  the  Landauer  equation  (2)  for  different  relative  unit-cell  sizes  (given  by  a).  This 
plot  shows  that  the  conductivity  plays  a  major  role  only  in  the  resonance  dumping  coefficient 
and,  in  general  for  small  conductivities  (as  expected  from  Born  approximation),  the  overall 
effect  on  the  magnetic  susceptibility  (x  in  the  lower  right  panel)  is  linear.  This  observation 
reaffirms  the  prediction,  based  on  the  scaling  principle,  that  the  resonance  frequency  of  the 
coil  mixture  is  mainly  a  function  of  the  geometrical  properties  of  the  unit-cell  while  being 
almost  independent  of  the  coil  conductivity. 

Figure  5  shows  the  extracted  resonance  model  parameters  as  a  function  of  the  number 
of  coil  turns  (two  upper  panels),  and  as  a  function  of  the  resonance  frequency  (two  lower 
panels).  The  middle  panel  shows  the  resonance  frequency  as  a  function  of  the  number  of 
coil  turns  and  can  be  seen  as  connecting  between  the  two  upper  panels  and  the  two  lower 
panels.  Three  curves  are  plotted  in  each  panel  for  three  different  relative  unit-cell  sizes. 
The  linear  dependence  of  the  curves  in  the  upper-right,  middle  and  lower-left  panels  as  the 
number  of  coil  turns  becomes  large  (or  resonance  frequency  becomes  small)  can  be  used  for 
simple  linear  extrapolation  (important  since  we  are  interested  in  the  THz  frequency  range 
-  below  10  THz)  that  is  shown  as  a  continuation  of  the  curves  as  dashed  lines.  The  linear 
extrapolation  of  the  upper-right  and  lower-left  panels  can  be  used  to  find  the  extrapolation 
of  the  curves  in  the  lower-right  and  upper- left,  respectively,  based  on  the  relation  between 
coil  turns  and  resonance  frequency  shown  in  the  middle  panel.  This  connection  between  the 
extrapolated  plots  in  different  panels  is  illustrated  in  the  figure  by  the  two  yellow  curved 
arrows,  beginning  with  the  panels  where  extrapolation  is  linear  and  crossing  through  the 
middle  panel  to  translate  turns  to  resonance  frequency  (and  vice-versa),  in  order  to  find  the 
extrapolation  at  the  panels  where  the  arrows  end.  The  linear  curves  for  a  large  number  of 
coil  turns  (or  small  resonance  frequency)  can  help  provide  understanding  of  the  nature  of 
the  relationship  between  those  values  and  the  resonance  properties. 

Figure  6  shows  the  peak-to-peak  variation  of  the  magnetic  susceptibility  (maximum  minus 
the  minimum  of  the  real  part)  normalized  for  the  surface  conductivity  of  one  conductance 
quantum  as  a  function  of  the  resonance  frequency  for  different  relative  unit-cell  sizes.  Such 
plots  are  meaningful  because  of  the  linear  effect  of  conductivity  shown  in  Fig.  4.  Again,  due 
to  the  linearity  of  the  curves  for  low  frequencies,  it  possible  to  extrapolate  them  down  to  the 
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FIG.  4.  Extracted  parameters  of  the  permeability  resonance  model  as  a  function  of  the  number 
of  modes  in  the  Landauer  equation  for  different  relative  unit-cell  sizes  (given  by  a);  (a)  resonance 
frequency,  (b)  dipole  strength,  (c)  damping  coefficient,  (d)  peak-to-peak  variance  in  magnetic 
susceptibility.  In  panel  (d) ,  the  dot-marked  curves  show  the  results  as  obtained  from  the  numerical 
simulation,  while  the  circle  markers  indicate  the  results  as  obtained  from  the  fitted  resonance 
model. 


THz  frequency  range.  This  plot  can  be  used  to  predict  the  magnetic  susceptibility  based  on 
a  desired  resonant  frequency  and  unit-cell  size,  for  a  given  surface  conductivity,  in  units  of 
conductance  quanta. 
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FIG.  5.  Dependence  of  the  extracted  resonance  model  parameters  on  the  number  of  turns  in  the 
coils  (two  upper  and  middle  panels),  and  on  the  resonance  frequency  (two  lower  panels). 


III.  LITHIUM  NIOBATE  (LN)  WITH  CIRCULARLY  ARRANGED  VACUUM 
HOLES 


Because  displacement  current  will  be  key  to  achieving  optical  magnetism,  we  consider  a 
structured  dielectric  material  that  should  provide  insight.  We  chose  LiNbC>3  -  LN  for  this 
initial  work  because  it  has  a  large  dielectric  constant  and  because  it  can  be  laser  machined 
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FIG.  6.  Peak-to-peak  variation  of  the  magnetic  susceptibility  normalized  for  the  surface  conduc¬ 
tivity  of  one  conductance  quantum  as  a  function  of  the  resonance  frequency  for  different  relative 
unit-cell  sizes.  The  dashed  curves  show  the  results  as  obtained  from  the  numerical  simulation, 
while  the  circle  markers  indicate  the  results  obtained  from  the  fitted  resonance  model. 

for  scaled  THz  measurements  and  imaged  using  light.  This  could  and  would  be  done  in 
collaboration  with  Prof.  Keith  Nelson’s  group  at  MIT. 

We  present  the  numerical  studies  of  extracted  homogenized  constitutive  parameters  er 
and  jir  from  structured  LN.  These  extracted  complex  parameters  are  given  as  functions  of 
frequency.  The  structure  under  study  consists  of  periodically  arranged  rings  of  vacuum  holes 
in  a  background  of  LN.  Each  unit  cell  measures  120  /irn  in  width  and  100  /irn  in  height.  In 
each  unit  cell,  seven  vacuum  holes  are  located  on  a  circle  of  30  /irn  radius,  equally  spaced. 
A  vacuum  hole  measures  5  /irn  in  radius.  The  structure  is  evident  in  the  electric  field  plot 
of  Fig.  7(a). 

We  investigated  the  homogenized  electromagnetic  response  of  the  structure  in  the  fre¬ 
quency  range  between  0.03  THz  and  0.3  THz.  At  0.194  THz,  the  wavelength  in  the  LN 
background  is  ALiNbo3  =  240  /irn,  twice  the  periodicity  of  the  crystal  structure,  which  results 
in  significant  reflection  (Fig.  7  b).  Frequency  regimes  meaningful  to  homogenization  should 
be  lower  than  the  first  crystal  mode,  and  in  order  for  the  effective  wavelength  to  be  greater 
than  10  times  the  structure  period,  the  frequency  ought  to  be  lower  than  0.039  THz. 
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FIG.  7.  (a)  The  in-plane  component  of  total  E- field  (incident  plus  scattered),  over  the  spatial 
extent  of  two  unit  cells.  A  TM-polarized  plane  wave  (H  out  of  the  page)  at  0.03  THz  is  incident 
from  the  left.  Each  unit  cell  is  composed  of  a  circular  array  of  holes  in  LN.  The  entire  structure 
under  study  consists  of  nine  such  unit  cells.  Shown  in  the  plot  is  a  region  over  two  unit  cells,  (b) 
Magnitudes  of  the  transmission  (S21)  and  reflection  (SR)  coefficients  for  an  array  of  holes  in  LN. 
The  first  crystal  mode  lies  at  around  0.2  THz.  The  structure  under  study  consists  of  9  unit  cells. 
Each  unit  cell  is  composed  of  a  circular  array  of  holes  in  LN. 


The  effective  constitutive  parameters  were  obtained  using  the  Nicolson- Ross- Weir  method 
[13],  and  are  given  in  Fig.  8.  Notice  that  for  frequencies  substantially  below  the  resonance 
(0.194  THz)  that  /ir  approaches  unity,  but  also  that  the  values  in  this  example  may  still  be 
useful  in  applications  (see  Fig.  8(d)).  The  point  of  this  example  is  the  develop  understand¬ 
ing  of  magnetism  based  on  displacement  current  (with  positive  real  part  of  the  dielectric 
constant)  that  could  be  translated  to  optical  frequencies.  This  result  suggests  further  ex¬ 
ploration  of  the  relation  between  the  geometry  and  homogenized  magnetic  response  should 
prove  interesting. 


IV.  COMPLEMENTARY  STRUCTURE  1:  CIRCULARLY  ARRANGED  LITHIUM 
NIOBATE  (LN)  RODS  IN  VACUUM  BACKGROUND 

The  homogenized  extraction  results  for  complementary  structures  having  dielectric  rods 
in  free  space  are  presented.  The  complementary  structure  to  that  used  in  Fig.  7  was  obtained 
by  replacing  the  LN  background  with  vacuum,  and  the  vacuum  holes  in  the  structure  by  LN 
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(a)  (b) 


(c) 
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FIG.  8.  Homogenized  results  for  holes  in  LN.  (a)  Extracted  relative  permittivity  er  over  the 
frequency  range  0.03  ~  0.3  THz.  (b)  er  at  frequencies  below  0.1  THz.  (c)  Extracted  relative 
permittivity  fir  over  the  frequency  range  0.03  ~  0.3  THz.  (d)  fir  at  frequencies  below  0.1  THz. 


rods.  Also,  the  radius  of  the  rods  was  increased  from  5  /irn  to  10  /irn.  The  reduced  distance 
between  rods  increases  the  coupling  and  enhances  the  ability  of  the  structure  to  scatter 
incoming  waves.  The  reduction  in  background  refractive  index  enlarges  the  frequency  range 
meaningful  for  homogenization.  Based  on  the  effective  wavelength  in  the  structure  being 
limited  to  be  no  less  than  10  times  the  period  (A  >  10 A),  we  estimate  that  meaningful 
homogenization  can  only  be  carried  out  at  frequencies  lower  than  0.25  THz.  Numerical 
scattering  parameter  magnitude  data  is  provided  in  Fig.  9.  The  extracted  homogenized 
responses  are  given  in  Fig.  10.  The  fine  features  in  the  extracted  effective  constitutive 
parameters  are  shown  in  the  magnified  plots  of  Fig.  11.  Significant  deviations  from  unity 
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have  been  achieved  in  the  homogenized  effective  magnetic  response  /ir ,  which  are  largely 
diamagnetic.  Note  that  the  deviation  of  fir  from  unity  is  larger  in  Fig.  11(b)  than  for  the 
complementary  structure  of  holes  in  LN  (Fig.  8(d)). 


FIG.  9.  Magnitude  of  the  S-Parameters  for  a  unit  cell  composed  of  a  circular  array  of  LN  rods, 
showing  the  first  crystal  mode  at  around  1  THz.  The  structure  consists  of  nine  unit  cells.  Each 
unit  cells  is  made  up  of  a  circular  arrangement  of  seven  LN  rods. 


Frequency  (THz)  Frequency  (THz) 


(a)  (b) 

FIG.  10.  Homogenized  response  for  the  complementary  structure  (LN  rods  in  free  space)  and  the 
data  in  Fig.  9.  (a)  Effective  relative  permittivity  er  over  the  frequency  range  0.03  ~  0.6  THz.  Both 
real  and  imaginary  parts  are  shown  in  the  plot,  (b)  Effective  relative  permeability  over  the 
frequency  range  0.03  ~  0.6  THz.  Both  real  and  imaginary  parts  are  shown  in  the  plot. 
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FIG.  11.  Expanded  scales  for  the  results  in  Fig.  10  for  LN  rods,  (a)  Real  part  of  the  effective 
relative  permittivity  er  over  the  frequency  range  0.03  ~  0.6  THz.  (b)  Real  part  of  the  effective 
relative  permeability  fj,r  over  the  frequency  range  0.03  ~  0.6  THz.  The  structure  under  study 
consists  of  nine  unit  cells.  Each  unit  cell  is  made  up  of  a  circular  arrangement  of  seven  LN  rods. 


V.  COMPLEMENTARY  STRUCTURE  2:  CIRCULARLY  ARRANGED  METAL¬ 
LIC  RODS  IN  A  VACUUM  BACKGROUND 

The  effect  of  replacing  dielectric  components  (LN)  with  metallic  elements  (having  negative 
real  part  of  the  dielectric  constant)  was  investigated.  The  structure  under  investigation  has 
the  same  unit  cell  and  periodicity  as  the  complementary  structure  in  the  previous  section.  In 
a  vacuum  background,  each  unit  cell  had  a  circular  arrangement  of  metal  rods.  Each  rod  was 
10  /irn  in  radius.  The  metal  was  a  hypothetical  material  having  permittivity  er  =  —41.7. 
A  field  solution  is  given  in  Fig.  12  and  the  calculated  scattering  parameter  magnitudes 
as  a  function  of  frequency  in  Fig.  13.  The  extracted  parameters,  er  and  /rr,  are  plotted 
as  functions  of  frequency  in  Fig.  14.  We  expect  the  homogenization  to  be  meaningful  at 
frequencies  below  0.25  THz.  The  effective  magnetic  response  is  mostly  diamagnetic,  having 
/ir  <  1.  Note  that  the  deviation  of  /j,r  from  unity  is  similar  in  Fig.  14(b),  where  er  <  0,  to  that 
for  the  material  having  the  same  geometry  but  with  er  >  0  in  Fig.  11(b).  The  relationship 
between  magnetism  and  the  sign  of  the  dielectric  constant  of  the  composing  material  needs 
to  be  investigated  further.  The  fact  that  we  achieve  similar  magnetism  independent  of  the 
sign  of  the  real  part  of  the  material  dielectric  constant  could  be  quite  important. 
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FIG.  12.  Metal  rods  in  free  space  (er  =  —41.7).  Plot  of  the  in-plane  component  of  E-field.  A 
TM-polarized  plane  wave  at  0.03  THz  is  incident  from  the  left.  At  this  frequency,  the  wavelength 
in  the  background  is  1  cm,  significantly  larger  than  the  size  of  the  unit  cell  120  /im.  Each  metal 
rod  has  radius  of  10  fim.  The  entire  structure  under  study  consists  of  nine  unit  cells. 


FIG.  13.  Amplitude  of  S-Parameters,  showing  the  first  crystal  mode  at  around  1  THz. 

VI.  A  COMPARISON  OF  DIELECTRIC  AND  METALLIC  STRUCTURES  UN¬ 
DER  TE-POLARIZED  INCIDENCE 

The  extracted  material  parameters,  er  and  /zr,  are  plotted  for  the  same  dielectric  and 
metallic  structures  that  were  covered  in  the  previous  two  sections,  but  with  TE  polariza¬ 
tion.  With  a  change  of  incident  plane  wave  polarization  to  TE,  we  were  able  to  probe  the 
other  tensor  components  of  the  constitutive  parameters,  allowing  investigation  of  differences 
between  the  two  complementary  structures  having  LN  rods  and  metal  rods,  although  un¬ 
der  TM  incidence,  the  effective  parameters  extracted  from  the  two  structures  were  almost 
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Frequency  (THz)  Frequency  (THz) 


(a)  (b) 

FIG.  14.  Homogenization  with  metal  rods  in  free  space  (er  =  —41.7)  for  TM.  (a)  Real  part  of  the 
effective  relative  permittivity  er  over  the  frequency  range  0.03  ~  0.6  THz.  (b)  Real  part  of  the 
effective  relative  permeability  fj,r  over  the  frequency  range  0.03  ~  0.6  THz.  The  structure  under 
study  consists  of  nine  unit  cells.  Each  unit  cells  is  made  up  of  a  circular  arrangement  of  seven 
metallic  rods. 

identical. 

In  Fig.  15,  the  S-parameters  of  these  two  structures  are  plotted.  The  incident  plane  wave 
has  TE  polarization,  with  the  electric  held  vector  oriented  parallel  to  the  rods.  Overall,  the 
metallic  structure  demonstrates  much  stronger  reflection,  and  heavy  attenuation  of  forward 
propagating  waves  than  does  the  LN  structure.  These  differences  in  scattering  characteristics 
lead  to  very  dissimilar  effective  medium  parameters,  as  shown  in  Fig.  16.  For  the  dielectric 
structure,  Fig.  16(a)  shows  er  >  1,  whereas  we  find  er  <  1  for  the  metallic  case,  as  shown  in 
Fig.  16(b).  The  results  agree  with  predictions  based  on  Maxwell- Garnett  mixing  theory. 

The  extracted  magnetic  responses  show  significant  deviation  from  unity,  as  plotted  in 
Fig.  17.  However,  we  need  to  consider  more  general  constitutive  relations.  While  spatial 
dispersions  can  give  rise  to  magnetism,  we  should  be  careful  in  reviewing  the  results  of 
Fig.  17  because  the  electric  held  (normal  to  the  page)  does  not  link  the  magnetic  hux,  except 
at  infinity.  The  extracted  response  should  in  fact  be  interpreted  in  terms  of  the  chirality 
tensor,  which  highlights  the  necessity  of  incorporating  material  chirality  or  bianisotropy  into 
the  extraction  methods.  We  note  that  the  inclusions  in  this  structure  do  not  exhibit  any 
handedness,  and  are  in  fact  non-chiral.  However,  such  a  lack  of  chirality  in  the  microstructure 
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FIG.  15.  (a)  Magnitude  of  S-parameters  for  the  complementary  structure  having  LN  rods. 

(b)  Magnitude  of  the  S-parameters  for  the  complementary  structure  having  metal  rods.  Both 
structures  consist  of  nine  unit  cells  with  the  same  cell  dimension. 

does  not  preclude  magneto-electric  coupling  and  could  provide  homogenized  parameters 
with  non-zero  entries  in  a  general  chirality  tensor  [17].  In  the  two  sections  that  follow,  we 
derive  aspects  of  plane  wave  propagation  in  isotropic  chiral  media  and  the  more  general 
bianisotropic  media,  in  anticipation  of  future  work  on  extraction  methods  that  takes  into 
account  these  general  constitutive  tensors. 


(a)  (b) 

FIG.  16.  Under  TE  incidence,  the  extracted  effective  permittivities  are  plotted  for  (a)  LN  structure 
(er  =  41.7  rods)  and  (b)  metallic  structure(er  =  —41.7  rods)  .  Both  structures  consists  of  nine  unit 
cells  with  the  same  cell  dimension. 
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(a)  (b) 


FIG.  17.  Extracted  effective  permeabilities  for  (a)  LN  structure  (er  =  41.7  rods)  and  (b)  metal¬ 
lic  structure  ( er  =  —41.7  rods)  .  Both  structures  consists  of  nine  unit  cells  with  the  same  cell 
dimension. 


VII.  PLANE  WAVE  PROPAGATION  IN  ISOTROPIC  CHIRAL  MEDIA 


Evaluation  of  magnetic  material  properties  requires  attention  to  the  material  symmetry, 
leading  to  the  constitutive  relations.  Given  these,  the  dispersion  relation  yields  insight 
into  wave  propagation  in  the  material.  This  picture  is  developed  here  and  will  be  used  in 
future  work  with  the  development  of  extraction  methods  incorporating  chirality  tensor  and 
anisotropy. 

The  plane  wave  eigenvalue  or  dispersion  surface  problem  in  relation  to  various  material 
tensors  (and  related  to  material  symmetry)  is  developed  here.  While  not  unique,  we  write 
D  =  /(E,  H)  and  B  =  /(E,  H).  From  Maxwell’s  equations,  this  has  a  fundamental  basis, 
because  J  drives  H,  for  example,  and  B  is  a  result  of  the  material  properties,  given  H.  This 
picture  is  quite  pervasive,  at  least  in  the  classical  field  literature.  However,  because  the 
homogenized  constitutive  parameter  relations  are  invertible,  the  specific  form  (D  =  /(E,  H) 
and  B  =  /( E,  H)  or  E  =  /(D,B)  and  H  =  /( D,  B))  is  largely  irrelevant,  it  would  seem. 
Therefore,  and  equivalently,  one  can  say  that  H  is  the  result  of  material  properties,  given 
B. 
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Consider  the  bi- anisotropic  constitutive  relations  of  the  form 


D  =  e0e  ■  E  —  jy/n 0e0  kt  •  H 

(4) 

B  =  jy/fi 0e0  R  ■  E  +  pop,  ■  H, 

(5) 

where  e  is  the  dielectric  constant  tensor,  fi  is  the  relative  permeability  tensor,  R  is  the 
chirality  parameter  tensor  relating  the  magnetic  held  to  the  electric  flux  and  vise  versa,  and 
exp (j cat)  time  dependence  has  been  assumed  [17].  The  form  of  (4)  and  (5)  implies  linear 
and  reciprocal  media,  and  is  convenient  in  the  solution  of  wave  propagation  problems.  This 
is  equivalent  to  what  has  been  called  the  Post  relations  (see  page  26  of  [17]).  The  crystal 
symmetry  dictates  the  form  of  the  tensors  in  (4)  and  (5)  -  see  Chapter  6  of  [17],  for  example. 

Removal  of  chirality  in  (4)  and  (5),  by  introducing  appropriate  cell  boundary  conditions 
of  through  appropriate  mixtures  leads  to  a  straight  forward  interpretation  of  e  and  Ji  from 
scattering  experiments  on  a  slab  of  this  material,  as  we  have  described  [12],  However, 
from  suitable  measurements  for  TE  and  TM  polarization,  all  the  constitutive  elements  are 
uniquely  available.  This  means  that  the  possibility  of  optical  magnetism  can  be  investigated 
through  optical  scattering  measurements. 

Consider  the  simplified  isotropic  and  chiral  constitutive  relations 

D  =  e0eE  +  j-H  (6) 

c 

B  =  — j-E  +  ^o/iH.  (7) 

c 

In  (6)  and  (7 ),  e  /i  and  k  are  unitless,  and  using  SI  units,  eo  has  units  F/m,  and  yU0  has 
units  H/m.  The  units  for  D  are  C/m2,  for  E  are  V/m,  for  H  are  A/m,  and  for  B  are  T  (or 
Wb/m2. 

For  a  plane  wave  with  spatial  variation  exp(— jk  •  r),  Maxwell’s  equations  become 


k  x  H  =  —  caD 

(8) 

k  x  E  =  caB 

(9) 

k-  D  =  0 

(10) 

k  B  =  0. 

(11) 

Note  that  the  reason  the  constitutive  relations  in  (4)  and  (5)  are  particularly  useful  is  that 
can  be  directly  substituted  into  the  right  hand  sides  of  (8)  and  (9).  A  forward  model  to 
describe  the  scattering  parameters  of  a  homogenized  slab  also  directly  use  this  form,  and 
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this  can  be  used  in  conjunction  with  measurements  to  extract  the  constitutive  parameters 
(and  evidence  of  magnetism)  [12]. 

For  a  plane  wave  with  exp(— jk  •  r)  spatial  dependency,  the  operation  Vx  is  equivalent 
to  —  jkx,  and  with  k  =  kxx  +  kyy  +  kzz,  we  can  write 

0  —  kz  ky 

k  x  =  |  kz  0  —  kx  ■  (12) 

0 


~ky  kx 


Using  (12),  we  can  write 


k  x  k  x  = 


p  _  p 

Ky  Kz 


kxky 


h  h  —  h2  —  P 

rvxrvy  rVy  rvz 


kx  kz 
kykz 


krkz 


h  h  _ b2  —  Jc2 

hvyhvz  t\jx  Ky 


(13) 


From  (9),  the  Helmholtz  wave  equation  is  formed  as 

k  x  k  x  E  =  cak  xB. 
and  substituting  (7)  and  then  (6)  leads  to 

k  x  k  x  E  =  cak  x  (  —  j  — E  +  ppH 


(14) 


c 


= —juj— k  x  E  +  ta/r0/u(— caD) 
c 

=  —  ju>— k  xE-wV  (e0eE+j-H 
c  V  c 

r,H  /  K 

=  -JU2 


K 


v  j-E  +  /Jo/iHj  -  cazpp  (^e0eE  +  j- HJ  (15) 

where  (9)  and  (7)  were  used  to  obtain  the  last  equation  in  (15).  Equation  (15)  leads  to 

k  x  k  x  E  +  fcg  (/re  +  k2)  E  +  j2ca2— /x0/iH  =  0.  (16) 

Starting  with  (8),  we  have 

k  x  k  x  H  =  —cak  x  D,  (17) 

and  substituting  (6)  and  then  (7)  leads  to 

k  x  k  x  H  =  —  cak  x(eoeE  +  j  — H 

V  c 

k. 


=  — caeoek  x  E  —  jca— k  x  H 
c 


K 


=  —  ca €q€  —  j  — E  +  pj/iH  )  —  jca— k  x  H 


=  — c a  eoe 


c 

-j-E  +  /J0/iH 

c 


■jca 


-j^K~,/u,/H)  .  (18) 
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Equation  (18)  leads  to 


k  x  k  xH  +  fcg  (jie  +  k2)  H  —  j2u2e0e—E  =  0. 


(19) 


Let  kx  —  0  and  ky  =  0,  i.e.,  consider  normal  incidence,  a  substantial  restriction  but  a 
simplifying  assumption.  Using  (16)  and  (19), 


A- 


Ex 

Ey 

Ht 


L  Hy 


=  0, 


(20) 


where  from  (13) 

— k2  +  /cg/re  +  k^K2 


A  = 


0 


0 


j2u2HQii\ 

0 


0 

j2u2fi0fj.* 


0 


— k2  +  fcg/ie  +  k^n2 
0 

— j2u2eo€^  0  — k 2  +  /cg/ue  +  k^n2  J 

We  modify  (20),  so  that  all  sub- matrices  have  the  same  units,  choosing 


-j2u;2e0e* 

0 


— k 2  +  /cg/re  +  k^K2 

0 


21) 


A'- 


Ex 
Ey 
Vo  Ha 

Vo  Hy  J 


=  0, 


(22) 


where  r) o  =  \J /i0 / eg  is  the  free  space  wave  impedance  with  units  of  and 


A'  = 


-k2  +  /cg(/ie  +  k2)  0  j'2K/ihg  0 

0  —  A;2  +  fcg(/xe  +  k2)  0  j'2/qufcg 

— j'2/qufcg  0  —  A;2  +  fcg(/xe  +  k2)  0 

0  —j2K/ihg  0  —  A;2  +  fcg(/xe  +  k2) 


23) 


All  entries  in  (23)  have  the  same  units.  The  dispersion  relation  (DR)  we  seek  is  detA'  =  0. 
We  can  use  the  identity 


det 


A  B 
-C  D 


=  detA  det(D  +  CA~lB) 
=  det  D  det  (A  +  BD^C) 


(24) 

(25) 
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or  other  analytical  methods  to  evaluate  the  determinant  of  (23)  and  find 

(det A')  ^  =  ^A'n  —  A/13A/31^ 

=  k4  —  2klk2(/xe  +  k2)2  +  k^(pe  —  k2)2 
=  0. 

Solving  (28)  as  a  quadratic  equation  gives 

k2  9 

~2  =  Ate  +  n2  ±  (2/cy //re) 
k0 

=  n 2, 

where  n  is  the  (complex)  refractive  index.  Consequently, 


(26) 

(27) 

(28) 


(29) 

(30) 


n  =  -yZ/ie  ±  /c.  (31) 

which  is  the  basis  of  chirality  leading  to  a  negative  refractive  index  [15].  However,  the 
specific  relationship  of  the  two  characteristic  solutions  to  the  polarization  must  be  considered. 
The  two  solutions  (n±)  have  been  assigned  to  right  and  left-handed  circularly  polarized 
eigenstates  [16],  attributed  to  [20].  A  treatment  in  terms  of  wave  fields  indicates  that  those 
are  circularly  polarized  and  have  the  refractive  indices  in  (31)  [21]. 

In  the  case  of  circular  polarization,  E  =  E0(x  +  jy )  for  RHCP  and  E  =  E0(x  —  jy )  for 
LHCP.  We  can  thus  write  Ey  =  pEx  with  p  =  j  for  RHCP  and  p  =  —j  for  LHCP,  for 
example.  When  |p|  ^  1,  the  polarization  state  is  in  general  elliptical  (linear  being  a  special 
case).  Consider  here  an  electric  field  basis  for  polarization,  i.e.,  use  of  E(f)).  However,  D 
has  been  used  [8],  where,  incidentally,  E  =  /( D,  B)  and  H  =  /(D,  B))  has  been  employed. 
The  description  of  D(t)  and  E(t)  differ  in  general. 

To  evaluate  the  polarization,  consider  (8)  and  (9)  with  the  constitutive  relations  in  (6) 
and  (7).  With  n  =  0  and  isotropic  material,  the  eigenvalues  for  k  are  identical  and  Ex  and 
Ey  are  independent.  However,  with  k/0,  the  situation  differs.  Writing  the  components  of 
(8)  and  (9),  with  k  =  kz,  and  rearranging,  gives 

kHy  =  uje0eEx  +  j~^Hx  (32) 

kHx  =  -ue0eEy  -  j—Hy  (33) 

c 

(jj 

hEy  j  Ex  uj  ft  ofiHx  (34) 

c 

(jJ  hb 

kEx  =  -j—Ey  +  upvpHy,  (35) 
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Arranging  (32)  -  (35)  into  matrix  form  yields  [7] 


-k 

— jiOhi/c 

Ex 

0 

capoR 

Hx 

—jum/c 

k 

Ey 

ta/iop 

0 

Hv_ 

-k 

— jcohi/c 

Hx 

0  ooe^e 

Ex 

■jcohi/c 

k 

Hy 

ue0e  0 

Ey 

Eliminating  the  magnetic  field  results  in 


(36) 


(37) 


Ex 

1 - 

O 

Ey 

1 

O 

(38) 


—j2k0kn  k 2  +  /cqK2  —  fcg/ie 
k 2  +  k^n2  —  /e,2pe  j2k0kn 

By  setting  the  determinant  of  the  coefficient  matrix  in  (38),  k/k0  —  n  —  y/pei/p  consistent 
with  (31).  Substituting  for  k  in  (38)  leads  to 


—A  =  ±j  for  n  =  y/JIe  ±  k,  (39) 

Ey 

i.e.,  yields  RHCP  and  —  n  corresponds  to  LHCP.  This  result  is  independent  of  the  sign 
of  k  and  holds  for  complex  /i,  e,  and  k. 

Consider  now  the  wave  field  transformations  [21]  that  allow 


D± 

=  ee0E  +  j-H  =  e±e0E± 
c 

(40) 

B± 

p/xoB  J  E  PitRo-B^,, 
c 

(41) 

which  are  provided  for  with  [21] 

E  =  E+  +  E_ 

(42) 

3V H  =  E+  -  E_ 

(43) 

H  =  H+  +  H_ 

(44) 

-E  =  H+  —  H_, 

<y~> 

(45) 

where  i]  =  Po\/r7c  with  p0  the  free  space  wave  impedance.  In  this  case, 


e±  =  ee0(l  ±  for) 

(46) 

/i-j-  flfl o(l  i  for)  i 

(47) 

where  nr  =  and  an  equivalent  non-chiral  medium  has  been  achieved.  Therefore, 

the  refractive  index  is  again  given  by  (31).  The  eigenmodes  E+,H+  and  E_,H_  have  been 
described  as  right  hand  and  left  hand  circularly  polarized  [21], 
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VIII.  PLANE  WAVES  IN  UNIAXIAL  BIANISOTROPIC  MEDIA 


In  preparation  for  developing  extraction  methods  that  incorporate  more  general  consti¬ 
tutive  tensors,  we  consider  now  the  description  of  plane  waves  in  media  that  have  both 
anisotropy  and  magneto-electric  coupling.  This  description  will  be  important  for  many 
nanostructured  materials  that  could  display  optical  magnetism. 

We  consider  reciprocal  bianisotropic  media  described  by  uniaxial  material  dyadics  [17]. 
Physically,  this  means  that  there  is  only  one  preferred  direction  in  space  defined  by  the 
geometry  and  the  positions  of  inclusions  or  molecules.  The  time  convention  taken  here  is 
exp(jut).  The  most  general  form  of  constitutive  relations  allowed  under  such  conditions  are 


D  =  e  •  E  -  j^eQHo(KtIt  +  «nzz  -  KJ )  •  H 

B  =  JJ  •  H  +  jyeo/ioMt  +  -  KJ)  •  E. 

where  the  permittivity  e  and  permeability  ~p  are  uniaxial  dyadics  given  by 

e  =  Co  (et/t  +  e«zz),  [1  =  +  Hn.zz) 


1  0  0 

0  0  0 

0-10 

0  1  0 

,  zz  = 

0  0  0 

,  j  =  Z  X  It  = 

1  0  0 

0  0  0 

0  0  1 

0  0  0 

The  cross  product  is  to  be  understood  as 


1 

<N 

CO 

1 

o 

bi 

a  x  b  = 

Q,  3  0  — Ol 

b2 

1 

o 

CM 

1 
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Therefore,  (48)  and  (49)  expand  into 


Dx 

1 

0 

0 

-M 

1  ^ 

Ex 

1 - 

0 

•+0 

£ 

Hx 

Dy 

= 

0  €-t  0 

Ey 

~  oho 

O 

1 

Hy 

Dz 

1 

e 

0 

0 

Ez 

- 1 

£ 

£ 

O 

O 

Hz 

Bx 

1 

0 

0 

Hx 

£ 

1 

O 

Ex 

By 

= 

0 

0 

Hy 

+  jVeoHo 

K  Kt  0 

Ey 

Bz 

0 

0 

Hz 

0  0  Kn 

Ez 

(48) 

(49) 


(50) 


(51) 


(52) 


(53) 


(54) 
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The  antisymmetric  elements  in  the  magnetoelectric  (chiral)  dyadic  arise  from  “hat”- 
inclusions  [17],  with  the  axes  oriented  along  the  2-axis.  In  the  formulation  of  Eqs.  (53) 
and  (54),  we  have  taken  the  optic  axis  of  the  uniaxial  medium  to  be  the  2-axis.  Note  that 
Ref.  [17],  page  28,  provides  a  proof  (to  be  verified)  that  Eqs.  (53)  and  (54)  satisfy  Lorentz 
reciprocity.  We  decompose  the  fields  as 

E  =  Et  +  En  z,  H  =  Ht  +  Hn  z  (55) 


We  take  the  two-dimensional  Fourier  transform  of  the  two  Maxwell  curl  equations  in  the 
transverse  (xy—)  plane,  so  that 


V  x  □  =  —  j 


0  dz  ky 

Wz  0  _  ^ 

hy  h‘x  0 


dD 

~dt 


=  jwD. 


(56) 


Maxwell’s  curl  equations  then  become 


-J 


—  —  k 

dz  ^ y 


£  o 


Ex 

En. 


=  -juj 


yt  0  0 

o  |U(  0 


Hx 

H„ 


+  K 


h' j  E  j- 
d't  Ey 


+  K 


Ey 
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Hy 
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K'tHy 
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Hx 

~ ky 

kx 

0 

Hz 

0 

1 

£ 

O 

Ez 

V 

i 

tu 

e 

£ 

0 

J 

(57) 


(58) 


where  k0  =  u>^/e0fi0  is  the  free-spaee  wavenumber.  Carrying  out  the  matrix-vector  multipli¬ 
cations  gives 


_9Ey_  I  .  p 

dz  '  ^y^z 

R tHx 

0 

0 

K't  Ex  E  Ey 

3 

EEa  —  k  e 

dz  ^X-^Z 

=  -ju 

0 

IkHy 

0 

+  k0 

KtEy  +  KEX 

kyEx  T  kXEy 

0 

0 

l-inH z 

3 

_dHJL  1  1  TT 

£t.Ex 

0 

0 
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dz  hjxIIz 
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0 

etEy 

0 

+ 

k0 

KHX 

kyHx  “j-  kXHy 

0 

0 

1 

Kl 

'sU 

- 1 

tq 

e 

(59) 


After  eliminating  the  normal  field  components,  Ez  and  Hz,  (59)  can  be  written  as  a  system 
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of  two  equations  in  the  transverse  field  components  as  (see  page  175  of  Ref.  [17]) 
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(61) 


Substituting  (60)  into  (61),  it  is  possible  to  eliminate  the  magnetic  field  and  yields  the 

lT 


second-order  wave  equation  for  the  transverse  electric  field  components  Et  = 
as 
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(62) 


where 
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with  n  =  en/in  —  k1 2 3  . 

Assuming  a  plane  wave  solution  with  the  longitudinal  (z-axis)  spatial  dependence 
exp(—  jj3z)  and  equating  the  determinant  of  (62)  to  zero  leads  to  the  dispersion  relation 
(pages  175  and  176  of  Ref.  [17]) 
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where  k2  =  k2  +  k2  and  nt  =  et/q  —  k2.  Equation  (64)  can  be  solved  for  /3,  giving  [17] 
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where  Re[\/D]  >  0  and 
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(66) 


After  solving  (65)  for  the  propagation  constant  /3,  the  cllipticities  (ratio  of  the  transverse 
held  components)  can  be  found  from  (62).  Note  that  in  finding  /3  =  (3'  +  j/3"  from  (65),  the 
sign  choices  need  to  be  carefully  investigated  based  on  causality  (Kramers-Kronig  relations) 
and  decay  of  the  Poynting  vector  for  a  passive  material. 
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